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R is shown that the dynamics of a single hole in a quantum antiferromagnet (described by the t-J 
model) can be simply understood in terms of a composite quasiparticle. This description provides 
naturally two different energy scales t and J corresponding to the inverse masses of the charge 
(holon) and spin (spinon) elementary excitations respectively. This picture is consistent with the 
exact results obtained on small clusters for the single hole spectral function and optical conductivity 
providing that one assumes the existence of a string-like force of magnitude J between the holon and 
the spinon. Then the hole quasiparticle can be interpreted as a bound state of its two constituents. 

PACS numbers: 74.72.-h, 71.27.+a, 71.55.-i 



I. INTRODUCTION 

The discovery of high-temperature superconductiv- 
ity has motivated a considerable effort in the study of 
strongly correlated fermion systems. The question of 
spin-charge separation is an important issue in these sys=. 
tems. In the case of the ID t-J model, it is knovifntl 
that a hole injected into the undoped system decays into 
two elementary excitations, a charge-1 spinless excita- 
tion (holon) and a neutral spin-1/2 excitation (spinon). 
Holons and spinous are independant excitations whose 
dispersion relations respectively scale with the hopping 
integral t and with the spin-spin coupling constant J, so 
that the decay products of an injected hole get quickly 
separated. The question of spin-charge seperation is 
more controversed in the case of the 2D t-J model, which 
is the simplest model used to desciibe the CuO planes 
of doped HTc materials. Andersoncl suggested that, in 
analogy with the 1-D case, a single hole in the 2D t-J 
model decays into elementaryj-excitations of spinon- and 
holon-type. Other approaches describe the single hole as 
surrounded by a region with reduced antiferromagnetic 
order, leading to a charge-1 spin-1/2 quasiparticle.- 

Based on exact calculations for small systemsQfl, we 
present here evidences that the hole in the 2D t-J model 
decays into elementary excitations of spinon and holon 
type. Following Ref. we make the case that the Drude 
conductivity and the total widthjjfJhe spectral density 
of the hole, which are observedElffl to depend mainly 
on parameter t, are inconsistent with the interpretation 
of the low energy behavior of the spectral density as a 
simple quasiparticle whose dispersion scales with J. All 
these properties of the hole can however be simply under- 
stood in terms of spinous and holons, provided that, in 
2D, there exists a long ranged attraction between spinons 
and holons which scales with J. In the presence of such an 
attraction, the quasiparticle-like pole seen in the spectral 
densityQU of the single hole can be regarded as a bound 
state of two constituents, a light one (holon) and a heavy 



one (spinon). 

The long ranged attraction between the heavy con- 
stituent and the light constituentjXtLthe hole in 2D can 
be visualized using a string picturetJ'tiJ. In this picture, a 
hole is initially created in an antiferromagnet by remov- 
ing an electron. The hole then hops away, leaving a wake 
of flipped spins behind. In absence of spin-flip terms in 
the Hamiltonian, the hole is bound to its initial position 
by the energy cost of the wake of flipped spins. A more 
formal description of this attraction can be obtained .iis,- 
ing the gauge theory of quantum antiferromagnetsllJEj, 
in which the fields of matter describing spinons and 
holons are coupled to a confining gauge fielcoEJ. This 
leads to the confinement of the constituents of the hole 
(spinon and holon), similar to the confinement of quarks 
which are the constituents of hadrons and mesons. The 
observation of spin-charge separation in the 2D t-J model 
at small time and length scales is therefore important, 
because it strongly suggests that elementary excitations 
with quantum numbers and interaction similar to those 
of quarks can be found in quantum antiferromagnets. 

In this note, we consider the properties of holes in the 
t-J model which is defined by the Hamiltonian 



<ji> 



(1) 



where the cinetic term T is given by 



T = Pc, 



-t 



<jl>a 



(2) 



and where Pq is the Gutzwiller projector which filters 
out states containing doubly occupied sites. The sum 
< jl > is performed over near-neighbor pairs, with each 
pair counted twice to maintain hermiticity. 
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II. SINGLE HOLE PROPERTIES 

Let us first consider the propagator of the hole, which 
is defined by 



1 



Cko 
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(3) 



where 



(TV) 



denotes the undoped (Neel) groundstate of 



the system with energy E^ ' , Ck 
and TV is the number of sites 

following, we set ?i = 1. In Fig. |^, we show the spectral 
density 



In Eq. (13) and in the 



Ak^ijj) = --IniGkaiuj) 

TT 



(4) 



obtained_.by exact diagonalization of 26-site and 32-site 
clustersO for J = and at different momenta. At 
E ((7r/2, 7r/2)), an 8t-wide continuum consistipg of two 
lobes with a hole in the middle is clearly visiblell3. When 
one goes along the F-M direction (from (0,0) to (tt, tt)) 
weight is re-distributed from negative to positive frequen- 
cies. This is consistent with the calculation of Ref. |l^ 
based on chiral spin liquid, in which these structures re- 
sult from the decay of the hole into a spinon-holon pair. 
As shown in Fig. |^(b,c), low energy poles located at the 
bottom of the spectrum become visible when the antifer- 
romagnetic coupling J/t is turned oil. They can be inter- 
preted as signatures of quasiparticleaJ with a definite dis- 
persion relation. Note that, although the quasiparticle- 
like pole is present for all momenta (only S is shown 
on Fig. 0) it has a much smaller weight at momentum 
k = (tt, tt) (M) as we shall discuss later on. 

Fig. ||(d-f) shows the ^sfoe quantities evaluated using 
a perturbative approachE3't2l based on linear spin-wave 
theory. In this approach, the fermionic operator Cka is 
expressed as a product of a spin-wave bosonic operator 
and a fermionic spinless hole operator. The hole prop- 
agator is then approximated by the propagator of the 
spinless hole, which is numerically pwaluated within the 
self-consistent Born approximationnjO'Eij. In this ap- 
proximation, the line of the spinless hole propagator is 
dressed by non-crossing spin-wave lines. The perturba- 
tive results of Fig. Hfd-f) agree well enough with the 
exact results of Fig.n2(a-c) to help to understand the 
structure seen in Fig. 0(a-c), as we shall discuss later on. 

The quasiparticle dispersion relation vs momentum is 
shown in Fig. ^(a) from an interpolation between the data 
obtained at the allowed K-values of the 26-cluster. It 
agrees remaxquably well with the Green Function Monte 
Carlo datao showing in particular a pronounced mini- 
mum at momentum (7r/2,7r/2) and a flat band in the 
vicinity of (7r,0). Fig. ^(b) shows the linear behavior of 
the total bandwidth defined as the difference in quasipar- 
ticle energy between the top of the band at F (expected 



to become exactly degenerate with M in the thermody- 
namic limit) and the bottom of the band at E. 
In Figs. 0(a-c), we show the weights 



(5) 



of the quasiparticle poles as a function of J/t for the 
bottom (Fig. |(a)) and the top of the band (Fig. |(b,c)). 
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where 



(JV-l) 



r J, ^ denotes the groundstate of the system 

with one hole and momentum fc. The vanishing of Z^. in 
the limit J/t ^ refiects the divergence of the "quasi- 
particle" size, i.e. the region of reduced antiferromagnetic 
order. 

At this point, it is important to notice that the spec- 
tral density extends on an energy interval of order 8 1 as 
seen in Fig. ^. It is remarkable that this energy width 
is basically almost independent of J while, on the other 
hand, the quasiparticle bandwidth varies linearly with J. 
As mentioned previously, the spectral function at higher 
energy also shows a strong k-dependence. These facts are 
clearly difficult to explain without taking into account the 
complex nature of the quasiparticle. 

The description of the hole propagator in terms of 
a simple quasiparticle dispersing like J is also a priori 
difficult to reconcile with the observation that the opti- 
cal conductivity of the hole depends mainly on t as first 
noticed by one of usQ. The optical conductivity is defined 
by 



N ^ En 



So 



(6) 



where the weight ttZ) of the Drude peak is simply pro- 
portional to the charge stiffness D given by 



D 
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where |0) and |n) respectively denote the groundstate and 
excited states with energies Eq and En of the system with 
a given number of electrons and where 



jx = it 



(8) 



with X denoting the vector connecting two neighboring 
sites along the x direction. 

Following Ref. |^ it is convenient to consider the ratio 
D/e^Uh since we expect D to scale with the hole doping 
n/i. This ratio has then the physical meaning of an in- 
verse mass. Fig. |^ shows the dependence of the quantity 
D / e^Uh given in units of i as a function of J/t for a sys- 
tem of 26 sites containing a single hole. Previous data 
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obtained by averaging over the boundary conditions (to 
reduce fisite size effects) of a 4 x 4 cluster with one or 
two holeg3 are also shown and are in excellent agreement 
with the -\/26 x -\/26 cluster data. The finite value ofpD 
reflects the ability of the system to conduct clectricityE^. 
Note the small dependence of D in J. It is difhcult to see 
how quasiparticles with properties that depend on J niay 
lead to a Drude conductivity which depends only on n3. 

III. SPINON-HOLON BOUND STATE 

We now make the case that the properties of the hole 
in the 2D t-J model can be simply explained in terms 
of elementary excitations similar to those of the ID t- 
J model. In the case of the 2D t-J model, spinons and 
holons can be formally defined within the gauge theory 
of quantum antiferromagnetscl as excitations dispersing 
respectively with J and t and which are bound together 
by a string potential V{v) — Q;J|r|, where a is a dimen- 
sionless constant. 

The mass of the charged excitation (holon) can be eval- 
uated from the response to an ptternal electromagnetic 
field by means of the f-sum rul£3 

/ duo(j^^{uj) ^ -K- — , (9) 
Jo 2m 

which relates the total weight of optical conductivity to 
the effective mass m and to the density n of charge car- 
riers. Using Eq. the quantity l/2m can be ex- 
pressed as a sum of two terms, (i) D/e^Uh shown in Fig. 
H (called "inverse optical mass" in Ref. |^) and (ii) the 
finite frequency integrated weight /q+ cfxx{'^) I (j^G^nh) 
also shown in Fig. H. This leads to 

m^^O.ri (10) 

in units where the lattice bound length is set equal to 
one, i.e. the holon mass depends mainly on parameter 
t. Note that this estimate for the holon mass agrees well 
with that of chiral spin liquid theoryH. Bgth Eqs. (||) 
and (|l^) are consistent with exact resultsEra for systems 
with several holes, in which the total weight is found to 
be nearly independent of J and proportional to doping. 

The quasiparticle-like pole appearing in the spectral 
density of Fig. ^ for finite J can be attributed to the 
spinon-holon bound state. The dispersion of the spinon- 
holon pair is characterized by the dispersion of the heav- 
iest object, which is the spinon in the case of interest 
J « {).2t. This is consistent with the finding that the 
quasiparticle bandwidth scales with J, shown in Fig. |^. 
Also, the vanishing of the quasiparticle weight Z in the 
limit J ^ seen in Fig. |^ reflects the divergence of the 
size of the spinon-holon bound state, consistent with the 
vanishing of the confining potential V{t) = a J|r| in this 
limit. 



As can be seen in Fig. ||(e,f ) , the spectral density eval- 
uated perturbatively for finite J shows a series of reso- 
nances above thcpttuasiparticle-like pole. These so-called 
string resonanceal3 can be interpreted as higher energy 
bound states of the spinon-holon pair. The energy de- 
pendence of both the quasiparticle-like pole and the first 
higher-energy resonance is displayed in Fig. ^. In this 
figure, the straight lines correspond to the fit 

En^pt + -i,,{J/tf/H, (11) 

where 13 = —3.28, 70 = 2.16 for the quasiparticle-like 
pole and 71 — 5.46 for the firS|tJughcr-energy resonance. 
This is the behavior expectedll^'O for the spectrum of a 
light particle in orbit around a heavy one, described by 
the Hamiltonian 

H = -^X/"^ + aJ\v\+ I3t, (12) 
2m 

where m is the mass of the light particle (holon) esti- 
mated using the f-sum rule. Using the energy depen- 
dence of either the quasiparticle-like pole or the first 
higher-energy resonance, one respectively findS|-a = 1.64 
or a = 2.16 for the string tension. Exact resultaHj for the 
quasiparticle-like pole are also well fitted using Eq.(|Tl|) 
with constants (3 = —3.359 and 70 = 2.77, leading to 
a — 2.39. Nevertheless the observation of high-energy 
resonances in exact diagonalizations of small clusters has 
been more controversial. Although some authorsQ have 
attributed some peaks in the spectral function of the 
4x4 cluster to these xesonances, similar studies carried 
out on larger clustersO for J > 0.3 failed to detect any 
sharp structure. However, it is tempting to attribute the 
structure appearing at a; ~ — 2.3i in Fig. ^(b) and at 
Lu « — 1.7t in Fig. |2|(c) (indicated by arrows on the plot) 
to the first higher-energy string resonance. This structure 
can be fitted using Eq.(|l|) with constants P = —3.359 
and 71 = 4.95, leading to a = 1.86. Note that these 
data correspond to quite small J/t ratios. For increasing 
J/t (let's say J/t > 0.2) the lifetimes of the string states 
becomes rapidly too small to be observable. 

As seen in Figs. |l| and || significant spectral weight ap- 
pears on a broad 8t-wide energy range. This continuum 
can be crudely viewed as the excitation spectrum of the 
light particle (holon) whose oscillator strength lies mainly 
at high frequencies. 

We conclude by a comment on the rather weak depen- 
dence of the Drude weight on J shown in Fig. |^ for one 
or several holes in small t-J clusters. This implies that 
a hole doped into the system has a dynamic mass which 
depends on t only. This is in contradiction with the lin- 
ear behavior of the quasiparticle bandwidth as a function 
of J seen in Fig. ||(b) , which implies that the mass of the 
hole depends on J only. It is very likely that one of 
these results is contaminated by finite size effects but, 
so far, it is impossible to determine which one without 
computations on larger samples. It is however important 
to settle the question of which of these results holds in 
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the thermodynamic limit in order to determine wether 
or not holes decay into more elementary constituents at 
low doping: If our observation that the Drude weight 
is independent of J is not seen on larger systems, this 
would mean that this decay is indeed not possible. This 
would be the case if the Drude "pole" separates, in the 
thermodynamic limit, into a legitimate Drude pole con- 
taining a small fraction cx J and a very low frequency 
band containing the rest. 

IV. ACKNOWLEDGMENTS 

DP wishes to thank the Stanford University for its hos- 
pitality during the early stage of this work and also ac- 
knowledges support from the EEC Human Capital and 
Mobility program under Grant No. CHRX-CT93-0332. 
DP and PB acknowledge IDRIS (Orsay, France) for al- 
location of CPU time on the C94 and C98 CRAY super- 
computers. 



^ W.F. Schmitt-Rink, CM. Varma and A.E. Ruckenstein, 
Phys. Rev. Lett. 60, (1988) 2793. 

C.L. Kane, P.A. Lee and N. Read, Phys. Rev. B39 (1989) 
6880. 

F. Marsiglio, A.E. Ruckenstein, S. Schmitt-Rink and CM. 
Varma, Phys. Rev. B43, (1991) 10882. 

G. Martinez and P. Horsch, Phys. Rev B44, (1991) 317. 
M. Boninsegni and E. Manousakis, Phys. Rev. B 46, 560 
(1992). 

23 W. Kohn, Phys. Rev. 133, (1964) A171. 
2^ R.B. Laughlin, preprint 

'^^ D. Poilblanc, H.J. Schuiz and T. Ziman, Phys. Rev. B46, 
(1992) 6435. 

FIG. 1. Single hole spectral functions for J=0 calculated 
on small 2D clusters with 32 and 26 lattice sites. The F, 
M and E points in reciprocal space corresponds to momenta 
(0,0), (7r,7r) and (7r/2,7r/2) respectively as shown in the in- 
sert. The spectral function at (7r/2, tt/2) is in fact obtained by 
averaging-the two spectral functions taken at the two nearest 
moment atj. 



^ P.A. Bares, G. Blatter and M. Ogata, Phys. Rev. B44, 

(1991) 130. 

2 P.W. Anderson, Phys. Rev. Lett. 64, (1990) 1839. 

3 J.R. Schrieffer, X.G. Wen and S.C Zhang, Phys. Rev. Lett. 
60, (1988) 944. 

* E. Dagotto, R. Joynt, A. Moreo, S. Bacci and E. Gagliano, 
Phys. Rev. B41, (1990) 9049. 

^ D. Poilblanc, H.J. Schuiz and T. Ziman, Phys. Rev. B47 
(1993) 3273. 

^ R.B. Laughlin, Conference Proceedings of the Summer 
School on Modem Perspectives in Many-body Theory, Can- 
berra, Australia (1993). 
D. Poilblanc, Phys. Rev. B44, (1991) 9562. 

* E. Dagotto, A. Moreo, F. Ortolani, D. Poilblanc and J. 
Riera, Phys. Rev. B45, (1992) 10741. 

^ D. Poilblanc, T. Ziman, H.J. Schuiz and E. Dagotto, Phys. 

Rev. B47, (1993) 14267. 
^° Z. Liu and E. Manousakis, Phys. Rev. B45, (1992) 2425. 
" E. Dagotto, Rev. Mod. Phys. 66, 763 (1994). 
^2 L.B. loffe and A.L Larkin, Phys. Rev. B39 (1989) 8988. 
" Z. Zou, J.L. Levy and R.B. Laughlin, Phys. Rev. B45, 

(1992) 993. 

M.C Diamantini, P. Sodano, E. Langman and G.W. Se- 

menoff, Nucl. Phys. B406, (1993) 1595. 

The diagonalization of the 32-site cluster is restricted to 

the most symmetric F and M k-points where use of the 

complete space group is possible. In this case the Hilbert 

space contains 37 596 701 symmetric states. 

On finite tilted clusters E is in fact slighly shifted from 

(7r/2,7r/2). This is a finite size effect d+ie to the constraint 

of tilted periodic boundary conditiona3. 

R.B. Laughlin, Phys. Rev. B45, (1992) 400. 



FIG. 2. Spectral function at for various J/t values cal- 
culated on the same 26-site cluster, (a-c) corresponds to the 
exact results and (d-f) to the spin-wave calculations. In (b,c) 
arrows indicate the location of the first string resonance. The 
dashed curves in panel (d-f) correspond to a 40 x 40-site clus- 
ter, for which finite size effects are negligible. 

FIG. 3. (a) Quasiparticle dispersion along some symme- 
try directions of the Brillouin zone obtained by interpolating 
exact results for a 26-sites cluster ( full line ) and obtained 
from spin-wave calculation for a 40 x 40-sites cluster (dashed 
line), (b) Quasiparticle bandwidth vs J/t calculated on a 
26-site cluster. Hexagons and triangles respectively denote 
exact and perturbative results. A linear behavior like 2.2 J/t 
is also shown for comparison. 

FIG. 4. Quasiparticle weights vs J/t calculated on a 26-site 
cluster for various location in k-space. (a), (b) and (c) corre- 
spond to the E, F and M points as indicated on the figures. 
Filled and empty squares respectively correspond to exact and 
perturbative results. 

FIG. 5. D/e^Uh for a single hole on a 26-site cluster (filled 
hexagons) in units oft. The finite frequency integrated weight 
/J^ o-xx{tjj) duj / {-KC^nh) is also shown (empty hexagons). Sim- 
ilar results are giv|«n for one (triangles) or two (squares) holes 
on a 4 X 4 clusteiu. In that case an average over the bound- 
ary conditions was used, leading to values larger than those 
obtained in Ref.B. 
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FIG. 6. Ground-state energy for a single hole as a function 
of J, obtained from exact calculations for a 26-site cluster 
(filled triangles) and from spin-wave calculations (empty tri- 
angles). Similar results are given for the energy of the first 
string resonnance (squares) at momentum corresponding to 
the band bottom of the quasiparticle-like pole. In the case of 
exact results, the first string resonnance is identified as the 
structure in the spectral density visible at u = —2.3t in Fig. 
||(b) and at Lj = -1.7t in Fig. ||(c). 
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